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Mixing in viscous fluids is challenging, but chaotic advection in principle allows 
efficient mixing. In the best possible scenario, the decay rate of the concentration 
profile of a passive scalar should be exponential in time. In practice, several authors 
have found that the no-slip boundary condition at the walls of a vessel can slow down 
mixing considerably, turning an exponential decay into a power law. This slowdown 
affects the whole mixing region, and not just the vicinity of the wall. The reason 
is that when the chaotic mixing region extends to the wall, a separatrix connects 
to it. The approach to the wall along that separatrix is polynomial in time and 
dominates the long-time decay. However, if the walls are moved or rotated, closed 
orbits appear, separated from the central mixing region by a hyperbolic fixed point 
with a homoclinic orbit. The long-time approach to the fixed point is exponential, 
so an overall exponential decay is recovered, albeit with a thin unmixed region near 
the wall. 

I. INTRODUCTION 

In many engineering applications, a viscous fluid must be blended with a substance, 
referred to as a passive scalar. This is generally called mixing, and it is well known that 
stirring greatly enhances this process. In fact, even if turbulence is unavailable (because of 
low Reynolds number or delicate substances), it is possible to mix rapidly by the process of 
chaotic advection [1, 2]. This involves the chaotic stretching of fluid particles by the flow, 
and the subsequent increase in concentration gradients of the substance to be mixed. The 
increased gradients facilitate the action of molecular diffusion, and homogeneity ensues. The 
net process (chaotic advection together with diffusion) is known as chaotic mixing. 

The quantity that is often tracked in mixing problems is the variance of the concentration 
of the scalar [5-12]. The variance is the spatial integral of the squared deviation from the 
mean concentration, and measures therefore the intensity of concentration fluctuations. The 
reasoning is that the variance tends to zero as the concentration is homogenized. Simple 
arguments [8, 12-19] suggest that the concentration variance should decay exponentially 
with time. In an idealized scenario the fluid particles are stretched exponentially and folded 
by the flow, yielding a characteristic filamentary structure as in Fig. 1. The filaments then 
achieves an equilibrium width where diffusion balances stretching [20]. Subsequently, the 
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FIG. 1. (a) Experiment with the 'figure-eight' stirring protocol, showing an advected blob of dye 
(India ink) in sugar syrup (details of the experimental apparatus are given in [3, 4]). The coordinate 
system used here is also indicated, as is the distance d(t) between the dark mixing pattern and the 
distinguished parabolic point on the wall. See Fig. 2 for the rod's path, (b) Numerical simulation 
for the same stirring protocol as in (a), but with the wall rotating at a rate Q = 0.4 rad/s. Here, d(t) 
is the distance from the mixing pattern to the hyperbolic point indicated by a dot. 

concentration field for that particle decays exponentially at the rate of stretching. An average 
over rates of stretching then gives the overall decay rate of the variance. In many cases the 
decay rate of the variance is determined in a less local manner [8, 12, 18, 21-24], but the 
decay is still exponential. 

This basic exponential-decay picture is appealing, but it is complicated by the presence 
of walls. In that case, several authors [25-34] have suggested that the no-slip boundary 
condition and the presence of separatrices on the walls slow down mixing: the decay rate is 
algebraic rather than exponential. Recent experiments [3, 4, 35] have confirmed this, and 
also showed that for a significant period of time the rate of decay of variance is dramatically 
reduced, even away from the walls, due to the entrainment of unmixed material into the 
central mixing region. 

In this paper, which is a more detailed and complete exposition of an earlier letter [36], 
we show that if the situation is such that mixing is slowed down by no-slip walls, then this 
can be cured by moving the wall in such a way as to destroy the separatrices. This creates 
closed orbits near the wall, effectively insulating the central mixing region from the wall (see 
Fig. 2). The rate of decay of the concentration variance becomes exponential rather than 
algebraic. The price to pay is that the thin region of closed orbits remains poorly mixed. Of 
course, in some applications moving or spinning the outer wall of a container is not practical, 
so we also discuss other mechanisms for creating closed orbits near the wall. 

The outline of the paper is as follows. In Section II we analyze the case of a device 
consisting of a stirring rod and a fixed outer wall. We describe the separatrices that appear 
at the boundary, and how the approach along the stable manifold of one of these separatrices 
is algebraic in time. In Section III we treat the same system, but with the wall rotating at 
a constant velocity. We find the wall separatrices are destroyed; instead, a hyperbolic fixed 
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FIG. 2. (a) Numerical Poincare section (stroboscopic map) of a single trajectory for the time- 
periodic figure-eight stirring protocol (dashed line). The circular domain has unit radius, the rod 
has radius 0.08, and the rod's two circular paths have radius 0.23. The chaotic region covers the 
entire domain, all the way to the wall, where two near-wall separatrices are indicated in the upper 
and lower portions of the domain, (b) Same parameters, but with the outer wall rotating at a rate 
of fl = 0.2 radians per period. Closed orbits appear near the wall, and the wall separatrices are 
destroyed. There is however a separatrix associated with a hyperbolic fixed point of the map (solid 
line), which isolates the chaotic region from the wall. The solid dot is the approximate position of 
the fixed point according to the model in Section III, and the cross is the numerically- measured 
fixed point. The discrepancy can easily be resolved by expanding to next order in f2. 

point appears away from the wall, with a homoclinic orbit separating the wall region from 
the central mixing region. The approach along the stable manifold of that fixed point is 
exponential in time. 

Section IV relates these results to the decay of concentration variance. We show that 
if the wall is rotating fast enough, the decay rate is not limited by the no-slip boundary 
condition, and mixing proceeds as rapidly as it would in the absence of a no-slip wall. We 
present numerical evidence for this in Section V, based on simulations of the evolution of the 
concentration field for a simple, but realistic, flow. We then offer some concluding remarks, 
and discuss other ways of creating closed orbits beyond physically rotating the wall. 



II. FIXED WALL 

Consider the experiment for the time-periodic 'figure-eight' rod-stirring protocol shown 
in Fig. 1(a). The numerically-computed Poincare section for the same 2-D Stokes flow 
(Fig. 2(a)) shows that this protocol promotes chaotic advection in the whole domain, as 
could be expected from the stretching and folding of dye filaments visible in Fig. 1(a). It 
also shows clearly that particles move very slowly near the wall, as evidenced by the closely- 
packed successive iterates, in accordance with the no-slip boundary condition. This suggests 
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FIG. 3. Solid line: The periodic function A{9) in Eq. (1) for the figure-eight stirring motion 
of Fig. 2(a). The two zeros correspond to the separatrices connected to the wall at the points 
highlighted in Fig. 2(a). The upper (stable) separatrix is at 9 = 9 S = tt/2, and the lower (unstable) 
separatrix is at 9 = 3tt/2. The dashed line is for the same stirring motion, but with the wall 
rotating at angular velocity f2 = 0.2 per stirring period. 



writing the flow near the wall simple map, 

9 = 9 + A(9)y + 0{y 2 ) , y = y- \A\9)y 2 + 0(y 3 ) . (1) 

where (6, y) is the position of a fluid particle at the beginning of a time interval, (9, y) 
its new position at the end. Here 9 is an angle measured counterclockwise along the wall, 
and < y <C 1 measures the distance from the wall (Fig. 1(a)). The first equation in (1) 
uses continuity and the no-slip boundary condition at the wall; then the second equation 
follows from incompressibility and the no-throughflow condition at the wall. The container 
is assumed to have unit radius. We can write (9,y) = (9(t),y(t)) and (9, y) = (9{t + T),y(t + 
T)), where t is time, T is the period, and t/T is the number of full periods of the stirring 
protocol that the system has undergone. If we truncate the map by neglecting the higher- 
order terms, then it preserves area only to first order in y. We shall see in Section III how 
to correct the map to ensure exact area preservation. 

Since it is continuous, the tangential velocity near the wall can only reverse sign at 
points 9 with A{9) = 0. These points correspond to the upper and lower wall separatri- 
ces (also called separation and reattachment points) visible in Fig. 2(a). (Note that these 
only act as separatrices near the wall, where the flow is almost steady; further away from 
the wall the manifolds of the separation points take a more convoluted form that permits 
chaotic transport in the bulk of the vessel.) There are only two such separatrices, so we 
conclude that A{6) must be as in Fig. 3, with the two zeros corresponding to the separa- 
trices in Fig. 2(a). (The functions in Fig. 3 were obtained numerically from 2D Stokes flow 
simulations [37].) 

Every point (9,y) = (#0,0) on the wall is a fixed point of the map (1). Letting (9, y) = 
(9o + 0, + Y), where (0, Y) are small expansion variables, we find the dynamics of particle 
trajectories near that point, 

9 = 6 + A(9 )Y + A'(9 )QY + O (<d 2 Y , Y 2 ) , (2a) 
Y = Y - \A\9 Q ) Y 2 + O (QY 2 , Y 3 ) . (2b) 
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The eigenvalues associated with the linearization of (2) are both unity, so (9q, 0) is a distin- 
guished parabolic fixed point [38] . Near most of these points the dynamics are uninteresting: 
particles just stream along the wall following 6 = + A(9 )Y, and approach or recede from 
the wall depending on the sign of A'(9 ). Eventually we must give up using (2), since the 
trajectories always leave the neighborhood of (6* ,0). However, for the two values of 6 for 
which A{9) vanishes (see Fig. 3), we get separatrices near the wall. 

We focus on the separatrix at 9 = 9 S = ir/2, where A(9 S ) = 0, and Eqs. (2) become 

= + A'(9 S )QY + O (0 2 F , Y 2 ) , (3a) 
Y = Y - \A\9 S )Y 2 + O (0F 2 , Y 3 ) , (3b) 

with (9, y) = (#s+0, 0+Y) and (0, Y) small expansion variables. Now the set {0 = 0, Y > 0} 
is invariant for small Y and corresponds to the separatrix, which is the stable manifold of 
the fixed point (9 S , 0). The evolution along the stable manifold is obtained by iterating 

Y = Y-\A\9 S )Y 2 1 (4) 

for small Y. This is a logistic map, and as we approach the fixed point at Y = we 
expect Y to change very little at each period. This suggests writing (4) as a differential 
equation for Y(t), 

^^^Y = -±A'(9 S )Y 2 , (F-F)/T«l. (5) 

The overdot denotes a time derivative. The solution is 

Y(t) = ^ (6) 

U l + lA'(9 s )Y(0)(t/T)- [) 

For this to represent the stable manifold, we require A'(9 S ) > 0, as is clear from Fig. 3. (The 
other separatrix exhibits finite-time escape to infinity, which takes particles away from the 
wall and into the bulk.) For long times, the rate of approach is 

Y(t) ~ (t/T)- 1 , t/T » {\A'{9 S )Y (0))" 1 . (7) 

The asymptotic form (7) for Y(t) is algebraic and independent of Y(0). The consequence of 
this independence is visible in the upper part of Fig. 1(a): material lines 'bunch-up' against 
each other faster than they approach the wall, thereby forgetting their initial position. 

The slow algebraic approach to the wall exhibited in (7) was shown in [3, 4, 35] to 
cause a drastic slowdown of the overall mixing rate in the vessel. Indeed, even though 
the fundamental action of the stirring protocol is chaotic, as evident in Fig. 2(a), a pool 
of unmixed fluid remains near the wall for very long times. More importantly, this pool 
contaminates the entire mixing pattern, all the way to its core, as unmixed fluid leaks 
slowly along the unstable manifold of the separation point at 9 — 37r/2, resulting in the 
characteristic cusped shape of the dye pattern in Fig. 1(a). This is because there is no 
barrier between the wall and the central mixing region. In the next section, we will see how 
we can create such a barrier by moving the outer wall. 
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FIG. 4. Iterates the area-preserving near-wall map map (9), which should be compared to the 
full numerical simulations of Fig. 2. (a) f2 = 0, with the upper stable distinguished parabolic 
fixed point shown as a dot; (b) f2 = 0.2, with the hyperbolic fixed point shown as a dot, and 
its associated homoclinic separatrix plotted as a solid line. We can clearly see closed trajectories 
outside the separatrix, near the wall. Once trajectories leave the vicinity of the wall, Eq. (9) no 
longer applies and we must look into the specific stirring mechanism. 

III. MOVING WALL 

Now we consider the case of a slowly-rotating wall, with mixing pattern as in Fig. 1(b). 
The near-wall map corresponding to (1) is 

6 = 6 + n + A(6)y, y = y -\K{Q)y\ (8) 

where < Q <C 1 is the angular displacement of the wall per period, and we have neglected 
higher-order terms in y. We assume that Q is small since in that case A{9) will change little 
from the case with Q = 0, as can be seen by the dashed line in Fig. 3. 

To help understand how the velocity Q modifies the wall map, it is instructive to plot 
iterates for a specific form of A(8). To match the numerical simulations, A(8) must be 
periodic in 6 and have only two zeros (Fig. 3); a simple model is then to take A{6) = — cos 9. 
For plotting iterates, it is preferable to use the exact area-preserving version of (8) (see 
Appendix A), 

e = e + tl + A(6)y } y = y -\A\6)y\ (9) 

which has unit Jacobian. (Using the exact map guarantees a faithful representation of closed 
orbits; but note that (6, y) are only approximate canonical variables, for small y.) Figure 4 
shows trajectories for the map (9). Once a particle leaves the vicinity of the wall, its 
trajectory becomes meaningless, since our expansions (1) and (8) are only valid for small y. 
However, the cusp structure for Q = in Fig. 4(a) is evident and is remarkably similar to 
the lower part of Fig. 1(a), where the unstable separatrix is located. 

Now we analyze the map (8) as we did in Section II when the wall was fixed. Again 
we look for fixed points of (8). All the distinguished parabolic fixed points on the wall 
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have disappeared, as well as the two separatrices. Since A(9) is continuous, has two zeros, 
and A'(9 S ) > 0, A(9) must have a minimum at 9i and a maximum at 62, where A' (6^2) = 
and hence y(0i,2,y) =y for all y. Seeking values for which the along-wall velocity also 
vanishes, we find there are fixed points at yi )2 = — f2/A(0i )2 ), which we require to be small 
for our analysis (that is, this defines how slow the rotation has to be). Since A(0 2 ) > 
(maximum) and A{9\) < (minimum), only 9\ has y\ > 0. (Recall that we assume Q > 0.) 
The other fixed point lies outside our domain. Hence, we focus on the unique near-wall fixed 
point —Q/A(9i)). (There could be others if A{9) has more extrema, but here we restrict 
to the case of two extrema.) 

We look at the linearized dynamics near the fixed point. Let (0, y) = (0i + 6, — fl/A(9 1 ) + 
Y); then 

= + A(9 1 )Y + 0(Q 2 ,Y 2 ,QY) , (10a) 
Y = Y - \A"{9 x )y\ + O(0 2 , Y 2 , GY) . (10b) 

To leading order in and Y, the motion near the fixed point is thus described by 

0\ _ / 1 A(9 1 )\ fe 



Y 



x -A»(0M l)\Y)- (U) 



The matrix in the above equation has eigenvalues 



»± = 1 ± ^- 1 2A(9 1 )A"(9 1 )y 1 = 1 ± y-^y n (12) 

where the argument in the square root is nonnegative since A(9i) < and A"(9i) > 0. 
For A"(9i) > and Q > 0, this is a hyperbolic fixed point, and the approach along its stable 
manifold obeys 

Y(t) ~Y(0) exp(-fit/T) (13) 
for (6(0), Y(0)) initially on the stable manifold, with decay rate 



to first order in Q. Compare this to (7): the approach to the fixed point is now exponential, 
at a rate proportional to the speed of rotation of the wall. In Fig. 5 we show the results 
for fx based on Eq. (14) and as measured in numerical simulations such as in Fig. 2(b). The 
two agree for small Q, as expected. In Section IV, we will see that the rate of exponential 
decay given by fx will dominate if it is slower than the mixing rate in the bulk. Otherwise, 
if fi is large enough, then the rate of mixing in the bulk dominates. 

Figure 2(b) shows a Poincare section for Q = 0.2. Very close to the wall (within a 
distance proportional to Q) trajectories are closed. A separatrix, consisting of a homoclinic 
orbit connecting the hyperbolic fixed point to itself, isolates the wall region from the bulk. 
Now it is the approach to this separatrix that will limit the decay rate, and as given by 
Eq. (13) this approach is exponential. Since the iterates of the map are close together on 
the separatrix, we can find an equation for the separatrix by writing a stream function 

V>(0, y) = \A{9)y 2 + ny , (u, v) = (dy^P , -cWO (15) 
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FIG. 5. The decay rate to the fixed point \x measured from simulations (solid line) and from 
Eq. (14). The two converge to each other for small Q. 

for a steady flow corresponding to the map. The value of the stream function at the hyper- 
bolic fixed point is tp(8i, —Q/A(6i)) = —^Q 2 /A(6i). The separatrix thus satisfies 



lA(9)y 2 s (6) + mO) = -^ 2 /A(e l ] 



which we can solve for y s (d), 



y s (9) = n 



1 + y/1 - A(6)/A(6 1 

MP) 



yi 



i - y/i - MpyMfh 

A(6)/A(6 1 ) 



(16) 



(17) 



where we have chosen the solution of the quadratic that ensures y s {@) > 0. The argument of 
the square root is always positive, since Q\ is a global minimum and A[Q\) < 0. When 9 = 9\ 
we recover y s {6i) = y±, as required. The singularity when A{9) = is removable, and we 
have then y s = \y\- The distance from the separatrix to the wall is proportional to Q; the 
point of closest approach to the wall is at 9 2 where A(8) is maximum, and the farthest from 
the wall is at Q\ where Aiff) is minimum (the hyperbolic fixed point itself). Figure 2(b) 
shows the separatrix and fixed point, computed using the theory in this section from the 
observed ^4(6*) for a moving wall (see Fig. 3, dashed line). The separatrix clearly divides the 
chaotic sea from the wall region. The cross in Fig. 2(b) is the numerically-computed fixed 
point for the map. The discrepancy between the numerical value and the theory is due to Q 
being relatively large in this figure (Q = 0.2), and can be easily accounted for by expanding 
to next order in Q. 



IV. DECAY OF VARIANCE 



So far in Sections II and III we have shown the following: (i) If the wall is not rotating, 
then there is a distinguished parabolic fixed point at the wall with a stable separatrix; the 
rate of approach of trajectories along that separatrix is algebraic, (ii) If the wall is rotating, 
the distinguished parabolic fixed point is destroyed, and instead a hyperbolic fixed point 
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appears elsewhere, away from the wall; the rate of approach of trajectories along the stable 
manifold of that fixed point is exponential. We shall now show that the different rates of 
approach are reflected in a different time-evolution for the scalar variance. 

The key is that the kidney-shaped 'mixing pattern,' by which we mean the shape of a 
blob that has been advected for several periods of the protocol approaches the fixed point 
at a rate d(t), where 

= i2T/(A'(9 s ) t), fixed wall [Eq. (7)]; 

U |d(0)e-^/ r , moving wall [Eq. (13)]. 1 ' 

(See Fig. 1 for the definition of d(t).) Indeed, a typical fluid particle first approaches the 
fixed point before being swept away along the wall or along the separatrix. 

Inside the central mixing region, we assume the action of the flow is that of a simple 
chaotic mixer. By this we mean that fluid elements are stretched, on average, at an expo- 
nential rate A. For simplicity, we also assume that in the absence of walls the concentration 
variance decays at the same 'natural' decay rate A, though in general the two rates can 
differ [15-17, 19]. Hence, after a time t a typical blob of initial size S will have length 5e xt . 
However, because of diffusion, its width will stabilize at an equilibrium between compression 
and diffusion at the Batchelor scale [11, 16, 17, 20] 



h = yfi/X, (19) 

where k is the molecular diffusivity. 

At every period, the pattern gets progressively closer to the wall. Assuming molecular 
diffusion can be neglected, area preservation implies that some white fluid must have entered 
the central mixing region. It does so in the form of white strips, clearly visible as layers 
inside the pattern of Fig. 1(a). In fact, if we assume that the mixing pattern grows uniformly, 
we can write the width A(t) of a strip injected at period t/T as 

A(t) = d(t) - d(t + T) ~ — T d(t) (20) 

where we also assumed that d(t) changes little at each period, consistent with experimental 
observations. This is positive since d is a decreasing function of time. 

Now, if a white strip is injected at time r < t, how long does it survive before it is wiped 
out by diffusion? The answer is the solution to the equation 

A(r) e- A (*- T > = £ B ■ (21) 

We interpret this formula as follows: The strip initially has width A(r) when it is injected; 
it is then compressed by the flow in the central mixing region by a factor e~ A ^~ T - ) depending 
on its age, t — r; and once it is compressed to the Batchelor length i B it quickly diffuses 
away. Thus, we can solve (21) to find the age of a strip when it gets wiped out by diffusion, 

t-r = \- 1 \og(A(r)/£ B ). (22) 

Eventually, at time t&, a newly- injected filament will have width equal to the Batchelor 
length. This occurs when 

A(t B ) = e B , (23) 
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which can be solved for £b given a form for A(t). After this time it makes no sense to speak 
of newly- injected filaments as 'white,' since they are already dominated by diffusion at their 
birth. Hence, the description we present here is valid only for times earlier than t-B, but 
large enough that the edge of the mixing pattern has reached the vicinity of the wall. 

In the experiment we measure the intensity of pixels in the central mixing region. We 
observe for 1 <C t/T ^ t&/T that the concentration variance is dominated by the amount of 
strips in the central region that are still white [3, 4]. Because of area conservation, the total 
area of injected white material that is still visible at time t is proportional to 

A w (t) = d(r(t)) - d(t) (24) 

where we use (22) to solve for r(t), the injection time of the oldest strip that is still white 
at time t. Hence, our goal is to estimate A w (t) for times 1 <C t/T < t-e/T, since A w is 
directly proportional to the concentration variance. To do this we need r(t), which requires 
specifying A(t). We examine the two possible forms in Eq. (18) in Sections IV A and IV B. 



A. Fixed Wall 

For the fixed wall of Section II, we have from (18) d(t) = 2T/(A'(9 s )t). The total area 
of remaining white strips at time t as given by (24) is proportional to 

■ 2T 2T 2T t-r 

From (18) and (20), the width of injected strips is A(t) = -Td = 2T 2 / (A 1 (9 s )t 2 ). Equa- 
tion (22) cannot be solved exactly, but since r(t) is algebraic the right-hand side of (22) is 
not large, implying that t/r ~ 1 for large t. We can thus replace r by t in (22) and the 
denominator of (25), and find 

• 4 -«^ M ^f M - SWT. (26) 

The decay of concentration variance is algebraic (~ l/t 2 ), with a logarithmic correction. 
The form (26) has been verified in experiments and using a simple map model [3, 4]. 



B. Moving Wall 

Now consider the case of a moving wall, with d(t) = d(0) e~ M * from Eq. (18). We 
have A(t) = -Td = /iTd(0) e"* 4 * = A(0)e"^. From (23), we have t B = /i" 1 log(A(0)/£ B ), 
and from (22), 

t-r = T ^(t B -t). (27) 

By assumption, r < t < t-&, so for consistency we require fj, < A, i.e., the rate of approach 
toward the hyperbolic fixed point is slower than the natural decay rate of the chaotic mixer. 
The area of white material in the mixing region is then obtained from (24), 

A w (t) = d(0) e""* (exp (y^B - *)) " l) , (28) 
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which in the regime X//i^$> 1 can be approximated by 

A w (t) ~ d(0) A~V (*b-*) e" M *, 



(29) 



The decay rate of the 'white' area is completely dominated by the rate of approach to the 
hyperbolic fixed point. The central mixing process is potentially more efficient (A > //), but 
it is starved by the boundaries. Hence, the wall slows down mixing, but unlike the fixed 
wall (Section IV A) the decay rate is still exponential. 

If the wall rotates rapidly enough to make /i > A, then the wall does not limit the 
decay rate at all. Indeed, for /i > A, we have t > t& in (27), since newly injected strips 
reach the Batchelor length £ B before strips that were injected previously. This violates our 
assumptions, and we conclude that in that case the white strips can be neglected; the decay 
rate of the concentration variance is then given by the natural decay rate A. 

In summary, we expect the variance decay rate to be given by 



that is, the decay rate of the variance is given by either the rate of approach to the hyperbolic 
fixed point //, or the 'natural' decay rate A, whichever is slowest. Because the decay rate 
is always limited by the natural decay rate A, one should not rotate the wall faster than is 
necessary to reach \x ~ A, as further increasing the rotation rate only decreases the size of 
the chaotic region, without any improvement to the mixing rate. We will see in the next 
section how this compares to simulations. 



The case of an algebraic decay was well-documented in [3, 4, 35], so we focus here on the 
exponential decay for a moving wall. Figure 6 shows the evolution of variance for several 
angular rotation rates Q. The variance was obtained by evolving a large number (of order 4 
million) of particles to model the concentration field of a passive scalar. This coarse method 
does not allow us to evolve the concentration field for a long time before particle statistics 
become inadequate, but we can conclude that the system enters an exponential regime 
for t/T > 25, with a rate of decay a. 

Figure 7 shows the numerically-measured decay rate \i to the hyperbolic fixed point near 
the wall (solid line), and the decay rate of the concentration variance (dashed line with error 
bars). For small Q, the exponential rate of decay is approximately equal to for larger Q it 
is less than //, since the decay rate is no longer dominated by the approach to the hyperbolic 
fixed point. In fact, as predicted by the theory, as the rotation rate is increased the decay 
rate tends to its "natural" value, approximated here by the dotted line. For small rotation 
rates, the agreement between \i and a cannot be called convincing, but is at least consistent. 
This is partly due to difficulties in measuring the decay rate accurately, as reflected by the 
large error bars on the dashed line. Note also that the effect is rather small here, so for this 
particular system the decay rate is not greatly limited by the walls. Nevertheless, we argue 
that the effect is there, and that in other situations this might be a more significant factor. 
For instance, it is known that the motion of walls can enhance heat transfer [39, 40]. 

Evidence that this type of effect could be important is provided by Figs. 8 and 9. They 
show the results of an experiment where the rod is moved along a "double-loop" rather 
than the figure-eight (see also [41] for a discussion of this system). The wall is fixed. Many 



a = min(/i , A), 



(30) 



V. 



SIMULATION RESULTS AND DISCUSSION 
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FIG. 6. Evolution of variance as a function of period, for different values of the rotation rate O. 
a is the fitted asymptotic decay rate (solid line). For small f2, a tends to fi — the decay rate 
onto the hyperbolic fixed point (see Eq. (14)). (The variance is measured over a finite subregion 
of the domain, which means it can increase at early times; the break in the top figure is due to an 
accidental absence of any particles in our measurement region, before particles get mixed.) 



industrial mixers, for example the wide class of planetary mixers [42, 43], impose such 
a multi-loop motion to paddles or stirrers. The double-loop motion, unlike the figure- 
eight, induces a net angular displacement of particles near the wall, leading to closed orbits 
(Fig. 8(b)). These closed orbits provide a similar effect to a moving wall: they prevent 
separatrices from being connected to the wall. Hence, the decay of concentration variance 
should be exponential, which is confirmed by the experimental results in Fig. 9. The analysis 
of such a system, however, is significantly harder than the moving wall case presented here. 
(Note that these two cases are not simply related by changing to a rotating frame: in general 
the rotation frequency of the wall and the period of the rods are not rationally related, so 
a fixed frame does not exist.) A detailed study of this case along the lines described here 
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(a) (b) 

FIG. 8. The 'epitrochoid' stirring protocol, (a) Experiment; (b) Poincare section, also showing the 
rod's trajectory. Notice the presence of closed orbits near the wall, even though the wall is fixed. 



would certainly shed some light on the transition from algebraic to exponential decay, and 
how system parameters could be tuned for improved mixing efficiency. 

Another important generalization is to extend the analysis to three dimensions, where 
even steady flows can exhibit exponential decay. This is particularly important for microflu- 
idics applications [44, 45], where the walls play an important role. The same analysis of 
separatrices at the walls should carry through, albeit with a richer range of possibilities that 
is inherent to three dimensions. 
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FIG. 9. Decay of concentration variance for the epitrochoid stirring protocol (Fig. 8), exhibiting 
an exponential decay. 
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Appendix A: Area-preserving Maps 

The near- wall maps (1) and (8) were derived by assuming a simple Taylor expansion 
in the vicinity of the wall, and applying boundary conditions. However, maps suffer from 
some drawbacks: they do not exactly preserve area, and they do not satisfy a natural time- 
reversibility formula. Though these drawback do not quantitatively modify the predictions 
of the paper, for completeness in this appendix we remedy the situtation. 

We wish to derive a symplectic (area-preserving) map for the flow having Hamiltonian 
(streamfunction) 

oo 

H(9,y,t) = uy + ±a(9)y 2 ]T 5(t-kT). (Al) 

k=— oo 

In this flow the fluid near the wall is rotated and 'kicked' at each period. The second term 
in (Al) is a perturbation of the integrable first term. Following Abdullaev [46], the simplest 
way to obtain an area-preserving map from (Al) is to introduce the generating function 

F{9 ) y) = {9 + n)y + \A{6)y\ (A2) 

with A(9) = a(9)T and Q = uT. With this type of generating function the mapping (9, y) — > 
(9, y) is defined implicitly via 

which gives the map 

9 = 9 + n + A(9)y, y = y-\A\9)y\ (A4) 
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Note that this map is the same as (8) to leading order in y. However, unlike (8) it ex- 
actly preserves area, by construction. One can see this directly from the generating func- 
tion F(9,y) = F(9,y(9,y)): from the second equation of (A3), we have 

QQ-V-Fee + Fey- ^ - - -— (A5) 



as well as 

?y = i = F d y = 1 

dy 6y dy dy F 9 y 

From the first equation of (A3), we have 



(A6) 



dd yo d9 yy d9 

as well as 



99 „ 09 dy F ee Fyy 

^^yy^-^ey ^ — , {*{) 



89 dy _ Fyy 

dy ~ dy ~ V (A8) 



Assembling all this, we find the linearization matrix 

8(9, y) 1 (Fh — FeeFyy F m 



d(9,y) F ey \ -Fee 



(A9) 



which always has unit determinant, as required for area preservation. The area-preserving 
map (A4) can be written in explicit form by solving for y in the second equation, 



y = Ul + 2A'{9)y - 1 /A' (9) , (A10) 



where the sign of the quadratic solution was chosen so that the map reduces to (8) for 
small y. 

The map (A4) solves one of the issues mentioned at the outset: it exactly preserves area. 
However, it fails to resolve the second issue, time-reversal symmetry. Indeed, interchang- 
ing (9, y) with (9, y) and letting T — > —T in (A4) leads to a different form for the map, 
which should not be the case for a map arising from a Hamiltonian such as (Al). To fix 
this is more complicated and requires the introduction of auxiliary variables (■&, T) [46]. We 
quote the form of the resulting map here for completeness: 

T = y - \A'{9)T 2 , $ = 9 + lA{9)T, (Alia) 
9 = $ + \A{9)T, y — T — \A'(9)T 2 , $ = $ + tt. (Allb) 

This map is both area-preserving and invariant under interchange of the barred variables 
together with T — > — T (the latter is equivalent to A — > —A and VI — > —Q). We can eliminate 
the T and d auxiliary variables in (All) by solving for T, 

T = 2 (yi + A'(6)y - l) I A! {9) , (A12) 

but the first equation in (Allb) must still be solved for 9, usually numerically. 



[1] H. Aref, 'Stirring by chaotic advection,' J. Fluid Mech. 143, 1-21 (1984). 



16 



[2] J. M. Ottino, The Kinematics of Mixing: Stretching, Chaos, and Transport (Cambridge Uni- 
versity Press, Cambridge, U.K., 1989). 

[3] E. Gouillart, N. Kuncio, O. Dauchot, B. Dubrulle, S. Roux, and J.-L. Thiffeault, 'Walls inhibit 
chaotic mixing,' Phys. Rev. Lett. 99, 114501 (2007). 

[4] E. Gouillart, O. Dauchot, B. Dubrulle, S. Roux, and J.-L. Thiffeault, 'Slow decay of con- 
centration variance due to no-slip walls in chaotic mixing,' Phys. Rev. E 78, 026211 (2008), 
http://arxiv.org/abs/0803.0709. 

[5] P. V. Danckwerts, 'The definition and measurement of some characteristics of mixtures,' Appl. 
Sci. Res. A A3, 279-296 (1952). 

[6] D. Rothstein, E. Henry, and J. P. Gollub, 'Persistent patterns in transient chaotic fluid mixing,' 
Nature 401 (6755), 770-772 (1999). 

[7] M.-C. Jullien, P. Castiglione, and P. Tabeling, 'Experimental observation of batchelor disper- 
sion of passive tracers,' Phys. Rev. Lett. 85, 3636 (2000). 

[8] D. R. Fereday, P. H. Haynes, A. Wonhas, and J. C. Vassilicos, 'Scalar variance decay in chaotic 
advection and Batchelor-regime turbulence,' Phys. Rev. E 65 (3), 035301 (R) (2002). 

[9] J. Sukhatme and R. T. Pierrehumbert, 'Decay of passive scalars under the action of single scale 
smooth velocity fields in bounded two-dimensional domains: From non-self-similar probability 
distribution functions to self-similar eigenmodes,' Phys. Rev. E 66, 056302 (2002). 
[10] G. A. Voth, T. C. Saint, G. Dobler, and J. P. Gollub, 'Mixing rates and symmetry breaking 

in two-dimensional chaotic flow,' Phys. Fluids 15 (9), 2560-2566 (2003). 
[11] E. Villermaux and J. Duplat, 'Mixing as an aggregation process,' Phys. Rev. Lett. 91, 184501 
(2003). 

[12] J.-L. Thiffeault and S. Childress, 'Chaotic mixing in a torus map,' Chaos 13 (2), 502-507 
(2003). 

[13] Y. B. Zeldovich, A. A. Ruzmaikin, S. A. Molchanov, and D. D. Sokoloff, 'Kinematic dynamo 

problem in a linear velocity field,' J. Fluid Mech. 144, 1-11 (1984). 
[14] B. I. Shraiman and E. D. Siggia, 'Lagrangian path integrals and fluctuations in random flow,' 

Phys. Rev. E 49 (4), 2912-2927 (1994). 
[15] T. M. Antonsen, Jr., Z. Fan, E. Ott, and E. Garcia-Lopez, 'The role of chaotic orbits in the 

determination of power spectra,' Phys. Fluids 8 (11), 3094-3104 (1996). 
[16] E. Balkovsky and A. Fouxon, 'Universal long-time properties of Lagrangian statistics in the 

Batchelor regime and their application to the passive scalar problem,' Phys. Rev. E 60 (4), 

4164-4174 (1999). 

[17] J.-L. Thiffeault, 'Scalar decay in chaotic mixing,' in J. B. Weiss and A. Provenzale, eds., 

Transport and Mixing in Geophysical Flows, Lecture Notes in Physics, volume 744, 3-35 

(Springer, Berlin, 2008), arXiv:nlin/0502011. 
[18] P. H. Haynes and J. Vanneste, 'What controls the decay of passive scalars in smooth flows?' 

Phys. Fluids 17, 097103 (2005). 
[19] P. Meunier and E. Villermaux, 'The diffusive strip method for scalar mixing in two dimensions,' 

J. Fluid Mech. 662, 134-172 (2010). 
[20] G. K. Batchelor, 'Small-scale variation of convected quantities like temperature in turbulent 

fluid: Part 1. General discussion and the case of small conductivity,' J. Fluid Mech. 5, 113-133 

(1959). 

[21] R. T. Pierrehumbert, 'Tracer microstructure in the large-eddy dominated regime,' Chaos 

Solitons Fractals 4 (6), 1091-1110 (1994). 
[22] A. Wonhas and J. C. Vassilicos, 'Mixing in fully chaotic flows,' Phys. Rev. E 66, 051205 



17 



(2002). 

[23] A. Pikovsky and O. Popovych, 'Persistent patterns in deterministic mixing flows,' Europhys. 

Lett. 61, 625-631 (2003). 
[24] D. R. Fereday and P. H. Haynes, 'Scalar decay in two-dimensional chaotic advection and 

Batchelor-regime turbulence,' Phys. Fluids 16 (12), 4359-4370 (2004). 
[25] S. W. Jones and W. R. Young, 'Shear dispersion and anomalous diffusion by chaotic advection,' 

J. Fluid Mech. 280, 149-172 (1994). 
[26] M. Chertkov and V. Lebedev, 'Decay of scalar turbulence revisited,' Phys. Rev. Lett. 90 (3), 

034501 (2003). 

[27] V. V. Lebedev and K. S. Turitsyn, 'Passive scalar evolution in peripheral regions,' Phys. Rev. 
E 69, 036301 (2004). 

[28] A. A. Schekochihin, P. H. Haynes, and S. C. Cowley, 'Diffusion of passive scalar in a finite-scale 

random flow,' Phys. Rev. E 70, 046304 (2004). 
[29] H. Salman and P. H. Haynes, 'A numerical study of passive scalar evolution in peripheral 

regions,' Phys. Fluids 19, 067101 (2007). 
[30] O. V. Popovych, A. Pikovsky, and B. Eckhardt, 'Abnormal mixing of passive scalars in chaotic 

flows,' Phys. Rev. E 75, 036308 (2007). 
[31] A. Chernykh and V. Lebedev, 'Passive scalar structures in peripheral regions of random flows,' 

JETP Lett. 87 (12), 682-686 (2008). 
[32] R. S. Mackay, 'A steady mixing flow with no-slip boundaries,' in C. Chandre, X. Leoncini, 

and G. Zaslavsky, eds., Chaos, Complexity, and Transport: Theory and Applications, 55-68 

(World Scientific, Singapore, 2008). 
[33] G. Boffetta, F. De Lillo, and A. Mazzino, 'Peripheral mixing of passive scalar at small Reynolds 

number,' J. Fluid Mech. 624, 151-158 (2009), arXiv:081 1.4519. 
[34] F. A. Zaggout and A. D. Gilbert, 'Passive scalar decay in chaotic flows with boundaries,' 

(2011), preprint. 

[35] E. Gouillart, O. Dauchot, J.-L. Thiffeault, and S. Roux, 'Open-flow mixing: Experimental 
evidence for strange eigenmodes,' Phys. Fluids 21 (2), 022603 (2009). 

[36] E. Gouillart, J.-L. Thiffeault, and O. Dauchot, 'Rotation shields chaotic mixing regions from 
no-slip walls,' Phys. Rev. Lett. 104 (20), 204502 (2010). 

[37] M. D. Finn, S. M. Cox, and H. M. Byrne, 'Topological chaos in inviscid and viscous mixers,' 
J. Fluid Mech. 493, 345-361 (2003). 

[38] L. E. Reichl, The Transition to Chaos: Conservative Classical Systems and Quantum Mani- 
festations, second edition (Springer, New York, 2004). 

[39] K. El Omari and Y. Le Guer, 'Numerical study of thermal chaotic mixing in a two rod rotating 
mixer,' Comput. Therm. Sci. 1 (1), 55-73 (2009). 

[40] K. El Omari and Y. Le Guer, 'Alternate rotating walls for thermal chaotic mixing,' Int. J. 
Heat Mass Transfer 53 (1-3), 123-134 (2010). 

[41] J.-L. Thiffeault, M. D. Finn, E. Gouillart, and T. Hall, 'Topology of chaotic mixing patterns,' 
Chaos 18, 033123 (2008), arXiv:0804.2520. 

[42] M. D. Finn and S. M. Cox, 'Stokes flow in a mixer with changing geometry' J. Eng. Math. 
41 (1), 75-99 (2001). 

[43] M. J. Clifford, S. M. Cox, and M. D. Finn, 'Reynolds number effects in a simple planetary 

mixer,' Chem. Eng. Sci. 59 (16), 3371-3379 (2004). 
[44] A. D. Stroock, S. K. W. Dertinger, A. Ajdari, I. Mezic, H. A. Stone, and G. M. Whitesides, 

'Chaotic mixer for microchannels,' Science 295, 647-651 (2002). 



18 



[45] E. Villermaux, A. D. Stroock, and H. A. Stone, 'Bridging kinematics and concentration content 

in a chaotic micromixer,' Phys. Rev. E 77 (1), 15301 (2008). 
[46] S. S. Abdullaev, Construction of Mappings for Hamiltonian Systems and Their Applications, 

Lecture Notes in Physics, volume 691 (Springer, Berlin, 2006). 



